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Abstract: We study the phenomenology of gauge singlet extensions of the Standard 
Model scalar sector and their implications for the electroweak phase transition. We de- 
termine the conditions on the scalar potential parameters that lead to a strong first order 
phase transition as needed to produce the observed baryon asymmetry of the universe. We 
analyze the constraints on the potential parameters derived from Higgs boson searches at 
LEP and electroweak precision observables. For models that satisfy these constraints and 
that produce a strong first order phase transition, we discuss the prospective signatures in 
future Higgs studies at the Large Hadron Collider and a Linear Collider. We argue that 
such studies will provide powerful probes of phase transition dynamics in models with an 
extended scalar sector. 



1. Introduction 



Explaining the origin of the baryonic matter of the universe remains an unsolved problem 
that lies at the interface of cosmology with particle and nuclear physics. Assuming that the 
universe was matter-antimatter symmetric at the end of inflation, the particle physics of 
the subsequently evolving universe would have to be responsible for the baryon asymmetry 
of the universe (BAU) observed in the present epoch: 



y = V± = J ( 7 - 95 ± °- 65 ) X 10_11 ' BBN i (j| 
B '~~ s 1 \ (9.29 ±0.34) x HT 11 , WMAP § 



where ng is the baryon number density, s 7 is the photon entropy density at freeze-out, the 
first value (BBN) is obtained from observed light element abundances and the predictions 
of standard Big Bang Nucleosynthesis, and the second value (WMAP) is extracted from 
the details of the acoustic peaks in the cosmic microwave background anisotropy power 
spectrum as probed by the WMAP collaboration. It is well-known that three ingredients 
must be present for successful baryogenesis ||: the violation of total baryon number (B); 
violation of both C- and CP- symmetry; and a departure from thermal equilibrium 1 . Al- 
though the Standard Model (SM), in principle, contains all three ingredients, the strength 
of CP-violation is too suppressed 2 and the departure from thermal equilibrium too gentle 
to bring about the observed value of Yg || . 

A variety of baryogenesis scenarios involving new physics have been considered over the 
years, with the energy scale of the corresponding dynamics ranging from the electroweak 
scale to the Planck scale (for recent reviews, see, e.g., || f?J). With the imminent start of 
operations at the CERN Large Hadron Collider (LHC), and the prospect for more powerful 
probes of new electroweak CP-violation using searches for the permanent electric dipole 
moments (EDM) of the electron, neutron, neutral atoms, and the deuteron || ^, 10 1 , it is 



a particularly interesting time to scrutinize the possibilities for electroweak baryogenesis 
(EWB), i.e. that the BAU was generated at the electroweak phase transition (EWPT) 
in the early universe. The EDM searches may tell us whether or not sufficiently large 
CP-violating asymmetries could have been created during the era of electroweak symmetry 
breaking, while LHC studies will provide information on the mechanism of EW symmetry 
breaking. In turn, this can lead to insights on the nature of the EWPT in the early universe, 
including whether, indeed, an EWPT occurred |ll| j. 

In this paper, we consider the implications of existing and future collider studies for the 
extensions of the scalar sector of the SM as they bear on the EWPT. In order to prevent 
the "washout" of any baryon asymmetry by electroweak sphalerons, the EWPT needs to 
be strongly first order - a requirement quantitatively characterized by the condition on the 
scalar (Higgs) background field cp = y/2(H°) 

^>1 , (1-2) 



The third ingredient can be avoided if CPT is violated. 
2 The strength of CP-violating effects in the electroweak sector of the SM is characterized by the Jarlskog 
invariant [II, J = (3.08±g;^) x 1(T 5 §. 
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where (H°) is the vacuum expectation value (vev) of the neutral component of the SU(2) l 
Higgs doublet and T c is the critical temperature at which the value of the minima of the 
Higgs potential at <p = and tp = ip c ^ become degenerate. Whether or not an EWPT 
satisfies Eq. (1.2) depends in detail on the scalar sector of the electroweak theory. In the 



SM, Eq. (|1.2p is equivalent to the following condition on the mass of the Higgs boson, m^: 



4E SM ( \ ) > 1 , (1.3) 



where vq ~ 246 GeV is the the value of <~p at T = 0, and TE$m is the calculable coefficient 
of the cubic term in the finite temperature SM scalar potential. For the value of E$m 
obtained in the SM, one finds that the Higgs mass must satisfy m/, < 42 GeV, in conflict 
with the direct search lower bound of 114.4 GeV from LEP 2 [fT^] . 

Extensions of the scalar sector of the SM can relax this requirement on . It is well 
known, for example, that in the Minimal Super symmetric Standard Model (MSSM), the 
presence of a light stop can enhance the cubic term, leading to Emssm ~ 10-EgM, allowing 



the mass of the lightest, SM-like Higgs to be as heavy as ~ 120 GeV [13], and thereby 



compatible with direct collider searches [ 14 ] . This enhancement arises when the soft mass 
Mjj is chosen to nearly cancel the finite temperature mass function n^(T). In this case, the 
order of magnitude increase in E arises largely from the number of stop degrees of freedom 
(= 2Nc)- Alternately, numerical studies have shown that extensions of the SM scalar 
sector containing singlet fields can also strengthen the EWPT [ TB|, It , [l?]] so that Eq. (|L 



is satisfied with scalar masses above the current direct search lower bound for the Higgs 



12]. Such singlets may be the low-energy remnants of scalar multiplets that are charged 
under gauge symmetries broken at high scales |l7|| . In supersymmetric models, they provide 
an attractive solution to the ^-problem, and models with dynamically generated /i-terms 
have been studied extensively in the literature [Jl8t| - A singlet field coupled to right-handed 
neutrinos can also give rise to a non-thermal production mechanism for sterile neutrino 
dark matter accommodating an explanation of both the cosmological dark matter and the 



observed velocities of pulsars [19] . 



In this paper, we analyze the generic implications of singlet extensions of the SM 
scalar sector for the EWPT and discuss the corresponding implications for low-energy phe- 
nomenology. Our goal is to identify generic features of singlet extensions that are relevant 
to the EWPT and to low-energy Higgs phenomenology, and that might arise in specific 
realizations, such as the above mentioned singlet extensions of the MSSM. To that end, 
we consider the minimal extension of the SM scalar sector containing one additional real, 



singlet scalar field, S [20|. Although this scenario does not contain additional CP-violation 
beyond that of the SM and, thus, does not directly address the new CP-violation needed 
for successful EWB, it does provide a framework for understanding general characteristics 
of the EWPT in models with extra singlet scalars. We will assume that the need for new 
CP-violation is addressed by aspects of a given SM extension that may or may not be 
directly related to the augmented scalar sector. 

In performing our study, we rely on a combination of analytic and numerical analysis. 
In doing so, we obtain an analog to Eq. (Ill|) that applies for phenomenologically viable 
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parameters: 

-e/T e 



V2 



cos a r 



2A 



c . 



+ •••>! , (1-4) 



where tana c is the ratio of the singlet vev and v/y2 at T c ; Te and e are coefficients, 
respectively, of the T-dependent and T-independent cubic terms in the finite temperature 
effective potential; A is the coefficient of the effective quartic term; Vq is the value of the 
degenerate minima of the potential at T c , defined with respect to the value at the origin; 
7 is a computable, positive constant that depends on the parameters in the potential; and 
the + • • • denote non-analytic or subdominant finite temperature contributions. 

From Eq. ( |1.4[ ), one sees that in the present context the strength of the EWPT can be 
enhanced with respect to the SM case in one of five ways: 

(a) increasing e compared to its SM value (esm = %V2Esm) 

(b) introducing new cubic terms in the potential giving a negative value for e 

(c) reducing the denominator through a tree level effects arising from quartic interactions 
of the singlet and SM scalar fields 

(d) reducing the denominator through singlet contributions to one-loop renormalization 
of the SU(2)l quartic coupling 

(e) generating a non-zero value for Vo just above the critical temperature. 

The first option generally requires introduction of a large number of additional scalar de- 
grees of freedom in order to be compatible with > 114.4 GeV, since the finite tempera- 
ture cubic term in the potential arises at loop level and is 1/4-7T suppressed. Similarly, option 
(d) requires either large singlet-doublet quartic interactions or multiple singlet scalars to 
overcome the one-loop suppression factors. Indeed, in the case we consider here of a single 
extra singlet degree of freedom, the singlet enhancements of e generally do not yield an 
appreciable strengthening of the EWPT (unless a large number of them are introduced, 



see Ref. [21| for a recent discussion of this case). In contrast, mechanisms (b) and (c) can 
yield a strong first order EWPT and phenomenologically viable scalar masses with inclu- 
sion of only one new singlet scalar having perturbative couplings to the doublet scalars, as 
they rely primarily on the presence of new tree-level terms in the potential. We also find 
numerous models that are compatible with low-energy phenomenology and that enhance 
the strength of the EWPT through mechanism (e). 

In addition to outlining the derivation of Eq. ( |l,4j ), we also discuss the low-energy 
probes of these mechanisms for enhancing the EWPT. In particular, we show that mech- 
anism (b) can be constrained by studies of Higgsstrahlung in e + e~ annihilation, since the 
tree-level cubic term introduces mixing between the SM and singlet scalars. We find that 
Higgs boson searches at LEP already impose stringent constraints on this mechanism, but 
for scalar masses lying above the LEP lower bound, there exists considerable room for 
inducing a strong first order EWPT in this way. Existing data from electroweak precision 
observables (EWPO) lead to additional constraints. It is well known that global analyses 
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of EWPO favor a rather light scalar 3 , and we find that models having a heavy singlet scalar 
that mixes strongly with the SU(2)-like scalar are disfavored. More generally, EWPO re- 
strict the mass of the SU(2)-like scalar to values below ~200 GeV while favoring small 
singlet masses and large mixing between light singlet-like and the SU(2)-like scalars. 

Future studies of Higgs boson production and decays at the LHC and a Linear Collider 
will provide powerful probes of scalar singlet models that generate a strong first order 
EWPT. For models in which the decay of one neutral scalar mass eigenstate into a pair 
of the second is kinematically allowed, one would expect deviations from SM Higgs decay 
branching ratios at the LHC. This possibility is particularly interesting in the case of 
mechanism (c) that is present in models having a Z2 (S <-> —S, H *-* —H) symmetry before 
spontaneous symmetry breaking (SSB). We find that the decay of the SM-like scalar into 
two singlet-like scalars is always kinematically allowed for models giving rise to a strong first 
order EWPT, and that the branching ratio to SM final states can be significantly reduced 
(from unity) for relatively light scalar masses. Thus, if a SM-like Higgs is discovered at 
the LHC, studies of its decay may yield important constraints on the viability of EWB in 
models with a Z2 symmetry. 

For models that do not have this symmetry, it is likely that either the singlet-like scalar 
can decay into two SU(2)-like scalars, or that neither scalar can decay into a pair of the 
other. One could probe the former class by searching for exotic final states at the LHC, 
such as four 6-jets or 6677. For models wherein no new scalar decay channels arise, future 
Higgsstrahlung studies at a 500 GeV e + e~ Linear Collider could provide an additional 
window on EWPT-viable scenarios. We discuss both possibilities below. 

Our analysis of the EWPT that leads to these phenomenological implications builds on 



extensive, previous work by others, both in the case of singlet scalar models [21, 24] and 



supersymmetric models with an extended scalar sector [|15| , |16| . In carrying out the present 
study, we have tried to amplify on earlier work by identifying general features of singlet 
scalar extensions that lead to a strong first order EWPT and that should apply to a broad 
array of specific models. We have also attempted to identify the ways in which electroweak 
precision measurements as well as present and future Higgs boson studies could be used 
to test EWPT-viable singlet extensions of the Standard Model scalar sector. Finally, we 
carry out a more detailed study of the pattern of symmetry-breaking associated with the 
singlet scalar, since models in which (S) / occur quite copiously. 

The discussion of our analysis is organized as follows: In Section |2| we outline the model 
for the singlet scalar extension of the SM, largely following the notation of Ref. |2(|. In 
Section [3| we analyze the pattern of symmetry breaking at temperatures for which SU(2) l 
is unbroken but the singlet scalar obtains a non-zero vev. Section ^ gives our analysis of 
electroweak symmetry breaking and the implications for the strength of the EWPT. In 
Section || we analyze the implications of low-energy (zero temperature) phenomenology 
for the EWPT in this scenario, including present and prospective constraints from e + e~ 



3 However, the forward-backward asymmetry for e + e — > Z° — > bb favors a heavy Higgs scalar, while 
the left-right asymmetries in e + e~ — > Z° — > hadrons indicate a Higgs with mass well below the LEP direct 



search bound Jz|. 
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annihilation, electroweak precision observables (EWPO), and Higgs production and decay 
at the LHC. We summarize our main conclusions in Section || 



2. Singlet Scalar Extension of the SM 



Following Ref. (2(J (and references therein), we study a minimal extension of the Standard 
Model (SM) scalar sector encompassing a single gauge singlet real scalar field S. The 
associated zero temperature, tree-level scalar potential reads 



where 



Vsm 
Vus 

V s = 



V = V SM + Vbs + V s 



2 H^h)+\ H^H 



= ~H 

°2 H2 , b 3 n 3 , °A 



(2-1) 



(2.2) 



S 2 



where H is the standard SU(2)/, scalar doublet with hypercharge one. Our notation fol- 
lows a practice sometimes used in the literature wherein no distinction is made between 
dimensionless and dimensionfull couplings (see, e.g., Ref. p0| ). Specifically, we use the 
symbol a for mixed singlet-SU(2)^ terms and b for singlet terms, and the subscripts refer 
to the corresponding power of the singlet field S. Tab. |I] should help the reader keep track 
of the various coefficients appearing in Vhs and in Vs, of their mass dimension and whether 
or not they are Z2 symmetric. 

We do not include a linear term in the singlet potential Vs, since for a generic potential 



V s = biS + 4$ 2 + + 



3 C 3 



2 3 

one can always remove the linear term by shifting S — > 
the equation 

h - b 2 f3 + b 3 f3 2 - b 4 p 3 



(2.3) 



5-/3, where (5 is a real solution to 
= (2.4) 



and obtain the 6j coefficients in Eq. (2.2) from 





b 2 = b 2 - 


- 26 3 /3 + 3b 4 p 2 






b 3 = b 3 - 


-3b 4 (3 






64 = 64 






Coefficient 


Corresp. Term 


Mass Dimension 


Z2 symmetric 


a\ 
a-2 
b 2 
b 3 
64 


(H^H) S/2 
(H*H) S 2 /2 

S 2 /2 

S 3 /3 

5 4 /4 


1 

2 
1 



No 
Yes 
Yes 
No 
Yes 



Table 1: A summary of the coefficients employed in the potential under consideration 
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In contrast to Ref. [20], we have not shifted S to remove its vev because we are 
interested in the pattern of symmetry breaking in the model at non- vanishing temperatures. 
Thus, we will keep the dependence on the neutral scalar vevs explicit. Letting and 
xq be the T = vacuum expectation values (vevs) of the neutral component of the SM 
Higgs H° and singlet scalar S, respectively, the potential minimization conditions 



dV 
dH 



(H°) T = =va/V2 



dV 
dS 



(2.5) 



{S) T =o=xo 



allow us to express two of the parameters appearing in Eq. (|2.2| ) in terms of the zero 
temperature vevs and the other parameters in the potential. For xq ^ 0, we choose to 
eliminate the mass parameters \j? and b 2 , 

fi 2 = X v 2 + (ai + a 2 x ) y (2.6) 

b 2 = -b 3 x -b,x 2 -^-^. (2.7) 

For xq = 0, the minimization condition with respect to S enforces the condition ai = 0, 
while fi 2 = Xov 2 as in the SM. In this case, as we discuss below, this is the necessary and 
sufficient condition for having a stable neutral scalar that can be the DM, as first noticed 
Ref. pq |. We emphasize, however, that imposing a tree-level 7j 2 symmetry on the potential 
(ai = = 63) does not imply a vanishing singlet vev. Only when xq = is it possible to 
have a stable neutral scalar. While this assumption is implicit in many previous analyses, 
we find that models with xq 7^ arise copiously in the present framework. 

The fields (h, s) describing fluctuations about the vevs are defined by H° = (v + h)/V2 
and S = xq + s, at T = 0. The corresponding entries in the mass matrix are given by 4 

d 2 V 

A = ^ = 2Ao^o 2 (2-8) 

/iJ^ = ^ + 2^-^ (2-9) 

2 d 2 V , . 

Vhs = ^r^~ = (°i + 2a2^o) v • (2.10) 
oaos 

The mass eigenstates h\ and h 2 are defined as 

h\ = sin 8 s + cos 9 h 

h 2 = cos 6 s - sinO h (2.11) 



where the mixing angle 9 is given by 

V_ 

1 + Vi + y 2 ' W_ pji-A*2 



Ian - V — = , where y = / fes g . (2.12) 



4 We discuss corrections resulting from the full Coleman-Weinberg effective potential below. These 
corrections lead to numerically small shifts to these conditions. 
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With this convention, | cos#| > l/y/2, therefore h\ is the mass eigenstate with the largest 
SU(2)-like component and h 2 that with the largest singlet component. The corresponding 
mass eigenvalues are given by 

m? )2 = ^^±^^v^T7 (2.13) 

where the upper (lower) sign corresponds to mi (777,2). 

For future reference it is useful to relate the parameters in V to those appearing in 
the notation of Ref. [ 20 1 , where the potential is written in terms of the zero-temperature, 
shifted field s only. One has 

V(H, s) = (h*h) + A {h^h) 2 + | (tftff) a (2.14) 

4W-(tf) s + ^ 2 + T s3 + T s4 ■ 

where 

61 = ai + 2a 2 x 
82 = a 2 

k 2 = b 2 + 26 3 x + 3b 4 xl (2.15) 

2 aiv 2 a 2 v 2 
= 6 3 x + 26 4 x - — - — 

k-z = b 3 + 36 4 x 
K4 = 64 , 

and where we have used the condition for a non-zero xo to eliminate 62 in terms of xo and 
the other parameters. From Eqs. (2.10- 2~l"5| ) we observe that even if the potential displays 



a Tj 2 (S «-> —S, H «-> — i?) symmetry before spontaneous symmetry-breaking, the zero 
temperature potential will not in general do so if the singlet vev xo 7^ 0. In this case, one 
would encounter mixing between the neutral SU(2)^ and singlet scalars, thereby allowing 
for decay of both mass eigenstates hi and h 2 to SM particles and precluding either from 
being a viable dark matter candidate. Only if 01, xq, and 63 all vanish does one allow for 
a (light) stable neutral scalar 5 . 

Various theoretical and phenomenological criteria restrict the values of the parameters 
in V: 

(i) Theoretically, V must be bounded from below. This condition requires the positivity 
of the quartic coefficients 64 and Ao- In addition, requiring the positivity of the quartic 
coefficient along an arbitrary direction implies that, if a 2 < 0, a| < 464A0. Alter- 
natively, one can require positivity of the determinant of the bilinear form including 
terms containing four powers of the fields, obtaining the same condition. 



'For non- vanishing 63, the singlet scalar can decay to a pair of SU(2)_l scalars at one-loop order. 



-7- 



(ii) Electroweak symmetry breaking is viable if (1) the determinant of the mass matrix 
is positive and (2) the electroweak vacuum is the absolute minimum of the potential. 
The first (necessary, but not sufficient) condition can be cast as 

b 3 x + 2b A x 2 -^ V 1 2 ' >0 , 2.16 

Ax 8Ao 

where v and x denote the vevs at arbitrary T <T C , with T c being the temperature at 
which the EWPT takes place. The second condition (that the electroweak minimum is 
the absolute minimum of the potential at T < T c ) cannot be cast in a simple analytic 
form: at a given temperature, the potential V in general has various minima. In our 
analysis below, we impose this condition numerically computing V at all minima and 
requiring that the electroweak vacuum is the true vacuum of the theory. 

(iii) Phenomenologically, searches for the Higgs boson constrain the combinations of pa- 
rameters that determine the masses mi2 and tan#. The LEP Higgs search results 
allow the existence of a scalar with mass below 114.4 GeV if its coupling to the 
Z-boson gHZZ is reduced from the SM coupling, g^zz- The corresponding Hig- 
gsstrahlung rate will be reduced by £ 2 = (gHZZ Iq^hzz) 2 compared to the SM rate if 
the decay branching ratios the same as in the SM. For the present case, this reduction 
factor is given by £f = cos 2 9 and £| = srn2 $ f° r the- two mass eigenstates h\ and /12, 
respectively. If the decay h\ — > /12/12 or /12 — ► h\h\ becomes kinematically allowed, 
then one has 

/_ \2 r SM 
c 2 _ ( 9HZZ \ 1 r91? s 

where T SM is the total decay width of of hi into conventional SM Higgs decay channels 
at the relevant mass. In our numerical scans discussed below, we will impose the LEP 
constraints in the (m^, £ 2 ) plane (k = 1,2), as given in Ref. |12[| . 

2.1 The finite temperature effective potential 

For purposes of analyzing the pattern of symmetry breaking involving both H and S, it 
is convenient to follow Pietroni |l5|] and work with a cylindrical coordinate representation 
of the classical fields v (T) and x(T). To that end, we define 93 = (p(T) and a = a(T) via 



v/V2 = <pcosa (2.18) 
x = (f sin a 

The resulting tree-level effective potential is 

V T= V, «) = (V cos 2 a - y sin 2 a) ip 2 (2.19) 

f a l 2 ^3 ■ 2 \ ■ 3 

+ I — cos a + — sin a I sin a tp 

, f\ 4 , a 2 2 -2 .^4.4^4 
+ I Aq cos a + — cos a sin a + — sin a j p , 
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where /j, 2 and 62 are implicit functions of the T = scalar vevs and of the other parameters 
in the potential as per Eqs. ( [2,6|JE7| ). 

The T = Coleman- Weinberg one-loop effective potential can be expressed in terms 
of the field-dependent masses rrii(v,x) = mi((p,a) pq]: 



G(y) 



^2n k G[m 2 k (ip,a)] 



(2.20) 



64tt 2 



1,1 ' W 



where Q is a renormalization scale (which we fix to Q = vo), the sum is over all fields that 
interact with the fields h and S, and n k is the number of degrees of freedom for the k-th. 
particle, with a minus sign for fermionic particles. The inclusion of V^ =Q leads to a shift 
in the minimization conditions 6 , corresponding to a shift in the dependence of /x 2 and 62 
on the zero temperature vevs and other parameters: 



1 d 



= ~ — V 1 

V ov 

1 d 



T=0 



(2.21) 



Abo 



dx ^ 



T=0 



Note also that the term proportional to —3/2y 2 in G(y) leads to a reduction in the mag- 
nitude of the quartic couplings Ao and 64 from their tree level values |23| . 

The finite temperature component of the effective potential receives two contributions: 
the one-loop component, and the bosonic ring contribution V r i ng that depends on 

the boson thermal masses. For the simple model we consider here, both contributions have 
been given in Refs. 27, 28] and we include those results here for completeness: 



2^ 2 ^ 

k 



^rikJB^F [ml(p,a)/T 2 ] 



(2.22) 



where the functions Jb,f{v) are given by 



Jb,f{v) 



dx x In 



1 =F exp ( — \J x 2 + y 



(2.23) 



withe the upper (lower) sign corresponding to bosonic (fermionic) contributions. Similarly, 
one has 

V, mg (ip,a,T) =-^-J2 n k {[M 2 (^a,T)f 2 - [m 2 (^a)f 2 } , (2.24) 



where the thermal masses M?((p,a,T) are given in terms of the T = masses and the 
finite temperature mass functions, 11^: 



M k (V, a, T) = ml(ip, a) + n fc 



(2.25) 



5 The correction to the mass parameters also occurs, as noted above. 
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Explicit expressions for the 11^ are given in Ref. J28|, and we do not reproduce them 
explicitly here. 

When T is large compared to a given mass m k , it is convenient to expand the functions 
Jb,f W.\ 



J B {y) * 


7T 4 

45 


vr 2 2 

H w - 

12 


n 3 

6 9 


2/ 4 
32 


J F (y) p 


. 77tA 
* 360 


"24 y " 


An 
32 


/V 

w 



In I — ] (2.26) 



where as = 167r 2 exp(3/2 — 2je), clf = vr 2 exp(3/2 — 2^e), and je is the Euler constant. 
The resulting expression for the full effective potential is given by 

V cS (ip, a, T) = Va(<p, a) + V?=°(tp, a) + AV(<p, a, T) (2.27) 

with 

AV(<p, a, T) = Vf # V, a, T) + V ling (tp, a, T) (2.28) 

90 720 y 

+ ^ Yl 9kml(ip,a) - — n k M%((p,a,T) 

k=B,F k=B 



—2 ^2 n k mt.(ip,a)]n 



64vr 2 

k=B,F 



a k T 2 



where g k = n k for bosonic degrees of freedom and g k = n k /2 for fermionic degrees of 
freedom and where Nb {Np) denote the total number of bosonic (fermionic) degrees of 
freedom, where the sum extends to the three Goldstone bosons as well. 

In the singlet extension of the SM of interest here, as in the SM, the transverse com- 
ponents of the W and Z bosons receive no finite temperature corrections to the masses, so 
that 

[M^,a,T)} T = ^ 2 = y cosW (2.29) 

rn#2/ rp\\T g 2 + g' 2 2 g 2 + g' 2 2 2 

[M z (ip,a,T)\ = v = cos a ip . 

In contrast, all other particle masses are screened by non-zero II- functions that are inde- 
pendent of the background fields and proportional to T 2 . Thus, the transverse compo- 
nents of the W and Z yield the dominant contributions to the effective cubic term in the 
Ves(ip, a, T) proportional to Tip 3 . In principle, it is possible to mitigate the screening of the 
scalar masses by a suitable choice of the tree-level parameters, thereby yielding additional 
contributions to the Tip 3 term. This strategy is similar to the case of the MSSM with a 
light stop, wherein the tree-level soft mass is chosen to cancel the thermal contribution. In 
the present case, however, we find that following a similar strategy for the thermal scalar 
masses generally leads to a potential unbounded from below, nonperturbative couplings in 
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ai/GeV 


0,2 


6 3 /GeV 


6 4 


xo/GeV 


A 


[-1000; 1000] 


[-i;i] 


[-1000; 1000] 


[0;i] 


[0; 1000] 


[0;i] 



Table 2: Ranges for the parameters used in our scans to generate the models employed for the 
figures in the remainder of the paper. 



the potential, or special, finely-tuned choices for the input parameters (such as xo). In a 
similar vein, Espinosa and Quiros (see Ref. p4|]), studied a related class of models contain- 
ing an additional complex scalar field, (j), having no vev and a potential that is invariant 
under <p ~~ * e ta (j). These authors observed that it is possible to choose the parameters in 
order to overcome plasma screening effects and enhance the strength of the EWPT. Nev- 
ertheless, they found that the maximum mass of the neutral SU(2) l scalar consistent with 
a strong first order EWPT is ~ 80 GeV when the dimensionless couplings in the potential 
are perturbative. Consequently, we will concentrate on the implications of terms in the 
T = potential that enhance the strength of the EWPT. 

We take here as the six free parameters the potential parameters a%, ai-, 63, 64, the 
quartic SM Higgs coupling Ao and the singlet vev at T = 0, xq. In the numerical scans we 
carry out in the present analysis, we linearly sample the above parameters in the ranges 
given m Tab. |. The ranges on the dimensionless parameters A, 64, and 02 are roughly 
consistent with the requirements of perturbativity as discussed in Refs. [24, 2S, pCf] . In the 
case of the Standard Model, these considerations imply that Xo/g 2 < 1, where g 2 « 0.4. 
One may obtain further restrictions on the value of 02 by requiring perturbativity up to 
a scale A that defines the limit of the validity of the theory, as has been implemented 
by Espinosa and Quiros in the case of the complex scalar model mentioned above p4fl . 
A detailed analysis of perturbativity considerations goes beyond the scope of the present 
study, where our focus falls on identifying general trends and for which limiting the ranges 
of the dimensionless parameters as in Table 1 should suffice 7 . 



3. The Singlet vacuum before the EWPT 

Recent studies of the pattern of EWPT in models with an extended scalar sector include 



both analyses where the singlet vev was assumed to be zero [£3, 21] and where the possibility 
of a non-singlet singlet vev prior to EWSB was studied numerically [28, 31]. In the present 
section, we address - both analytically and numerically - the role of a non- vanishing singlet 
vev on the strength of the EWPT. 

Spontaneous symmetry breaking in the singlet sector offers two possible scenarios: (i) 
the singlet scalar has already acquired a vev before the EWPT (i.e. at T > T c ), or (ii) 
the singlet has zero vev. We show below that under reasonably general assumptions, it 
is possible to derive an analytical condition on the parameters appearing in the potential 



7 As we also discuss below, imposing constraints from EWPO restricts the scalar masses to be rather light, 
corresponding to relatively small magnitudes for the dimensionless couplings in the potential. Choosing 
these couplings near the upper ends of the ranges in Table 1 generally leads to scalar masses that are 
inconsistent with EWPO. 
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to determine which of the two options occurs. This is necessary in order to determine 
the value of the EWPT critical temperature, T c , that is defined as the temperature at 
which the minima of broken and unbroken electroweak symmetry become degenerate. It 
is possible that the minimum of unbroken electroweak symmetry at T > T c corresponds 
to a non-vanishing singlet vev, so it is useful to determine the conditions under which this 
possibility occurs as well as the value of the corresponding unbroken minimum. 

To that end, we consider the effective potential for a = tt/2 for which the SU(2)l vev 
is zero. It has the general form 

y eff (y?,7r/2,T) = 6 + 61 99 + y<^ 2 + y <^ 3 + j^ 4 + • • • , (3.1) 
where the bi are T-dependent functions of the tree-level parameters 

6 = coT 4 + doT 2 (3.2) 

— ! :•■ I 1 - -±-± ) 01 + (l - 01 (3.3) 




b 2 = h + ^tSl 1-^-^ ) a 2 + ( l-h/^± ) a2 (3.4) 



7T 



63 = 63 (3.5) 

64 = 64 ; (3-6) 

where the "+ • • •" denote non-analytic contributions arising from the logarithmic field de- 
pendence in V\\ where the dependence of cq and a*o on the parameters in the potential is 
not important for this discussion and, therefore, not shown explicitly; and where we have 
not explicitly included the analytic tp 2 and (p 4 contributions arising from V^ =0 that are 
suppressed by 1/64-7T 2 . The quantities yfai have been obtained by expanding the thermal 
masses Mk(<p, tt/2, T) in powers of mk(<p, n/2)/T to second order and are given by 

a x = 9 - + ~^ + y l + ^ (3.7) 
x 4 2_ 4 12 y ' 

_ 9 2 3A yf a 2 

a+ 4 2 4 12 

h , a 2 

ot- = 

4 3 

The potential in ( |3.l|) can in principle develop a global minimum at cp 7^ 0. In order to 
analyze the latter possibility, we first make the assumptions that (a) the non-analytic in- 
dependence arising from the screened masses and logarithmic term in V 1 T=0 have negligible 
impact, and (b) the potential has a Z 2 symmetry or that the breaking of this symmetry is 
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sufficiently weak that the b\ip term in Eq. (|3.1| ) is negligible (for analytic ease we shall set it 
to zero in what follows). Under these assumptions, the equation dV e s((p, a = i:/2)/dtp = 
has the following solutions: 

<p = (3.8) 

*± = \ ^ 

where the last solution is physical as long as A = 63 — 46264 > 0. In what follows, we 
discuss the algebraic conditions that dictate which is the global minimum of the potential 
before EW symmetry breaking, i.e. whether or not tp ^ at T > T c . 

If A < the origin is the high T minimum, since d 2 V e g((p, a = n/2)/dp 2 = 62, 
62 > 63/(464) and 64 > 0. If A > the origin is a maximum if 62 < (case 1.), and a 
minimum if 62 > (case 2.). 

1. If 62 < 0, <p± are the two minima of the potential, and 

/ 2 >-(|+N (!-¥)• <«»> 

To establish which is the true minimum, we compute 

VA (b\ 



V cS {ip + ) - V eS {y-) = |p- - 6 2 J 63 (3.11) 

from which one reads (recalling that since A > and 64 > the term in parenthesis 
is positive) that if 63 > the minimum is (f-, while if 63 < it is </?+. 

2. If 62 > 0, the minimum is either the origin or ip + (</?-) for 63 < ( 63 > 0). The 
origin is the absolute minimum if V eS (ip±) > i.e. (substituting (|J) into ( |3TTo| )) if 

66264 - b\ > =f6 3 \/A for <p± (3.12) 

From the equation above, simply dividing right and left hand sides by 63, it follows 
that the origin is the absolute minimum if 

6 2 6 4 2 

if Slr >- . (3.13) 

In contrast, for R < 2/9 and 63 > the high-T minimum is at <f—, while if R < 2/9 
and 63 < it is ip+. 

In summary, in the limit where 6q = b\ = 0, the global minimum of the potential 
before EW symmetry breaking lies at the origin if R > 2/9 and at ip 7^ if R < 2/9. 

When discussing the strength of the EWPT, it will be useful to know the magnitude 
of the high-T minimum before EWSB. It is straightforward to show that for < it! < 2/9 

V min (<p min , a = it/2, T) = Vo = -^pf(R) (3-14) 

f(R) = - (1 - 4R)(1 - 2R) - 2R 2 + (1 - 4i?) 3/2 
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The function f(R) decreases monotonically from one to zero as R is increased from zero 
to 2/9. For R > 2/9, Vq = 0. Thus, if 63 is of order the electroweak scale or smaller 
and 64 ~ O(l), then one would expect from Eq. (|3.14j ) that |Vb|/T 4 << 1 for T of order 



the electroweak scale. As we discuss below, our numerical study suggests that the critical 
temperature for the EWPT, T c is typically ~ 100 GeV for phenomenologically viable choices 
of the potential parameters. Thus, we will adopt an expansion in | Vo | / T^f in analyzing the 
strength of the EWPT. Our result in Eq. ( |1.4[ ) relies on the validity of this expansion, that 
we have confirmed numerically. 

The foregoing analysis can be generalized to the case where the linear term in V e s((p, vr/2, T) 
is not negligible by shifting the classical field 99 by a T-dependent constant (3{T) as in 
Eq. ( |2.4D to eliminate the 61 term and defining a new set of coefficients as in Eqs. (]2.5|), 
The conditions for a nontrivial high-T singlet vev that has Vo < will be the same as 
above, but expressed in terms of the shifted coefficients. The connection with the param- 
eters in the tree-level potential will be implicit in these relations, which can be studied 
numerically. We have carried out such a numerical investigation of this issue, and we find 
that for models which do not exhibit a TL 2 symmetry and that are compatible with LEP 
bounds on (mk,Ck) (k = 1,2), one has |Vo|/T 4 <C 1. 

We checked the validity of our criterion outlined above against the full numerically 
computed finite temperature potential. We find that for almost all models with R > 2/9 
the origin is the absolute minimum at T > T c ; however, we do find a few models with 
< R < 2/9 where the origin is the absolute minimum, indicating that R < 2/9 is a 
necessary but not sufficient condition for the singlet field to develop a non- vanishing vev 
at T >T C . We also numerically compared the quantity R in Eq. ( p. 13 ) with the quantity 



Ro = 62&4/&I defined with the tree level couplings appearing in the potential, assuming a 
vanishing linear term b\ in the potential, and we find that the deviation of Rq from R at 
T = T c is typically smaller than a few percent. Henceforth, for practical purposes, Rq can 
be used to discriminate potentials where the singlet acquires a vev before the EWPT from 
those where it does not. 

4. The electroweak phase transition 

In order to explore the implications of the augmented scalar sector for the EWPT, it is 
convenient to express the effective potential for general values of a and (p as 

Ves(<P, a, T ) = V + (e - eT)ip 3 + [2D(T 2 - T 2 ) + (b 2 sin 2 a)/2] (p 2 + BT 2 tp-\ (4.1) 

where 

f ai + b 3 \ . 

B = sin a 

V 12 J 

D = Dsm cos 2 a + — [02(1 + sin 2 a) + 364 sin 2 a] 
48 

e = 2V2E SM cos 3 a 

( ai 2 , ^3 • 2 A ■ /i n\ 

e = I — cos a + — sin a I sin a (4.2) 
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A = Asm CO cos 4 a + ^ cos 2 a sin 2 a + -^ sin 4 a + 



where -Dsm, ^SM and Asm correspond to the usual SM values [3C]; and where the "+ • • •" 
denote terms involving non-analytic dependence on ip associated with the one-loop finite 
T-effects and the screened thermal masses. In principle, one could expand the latter in 
powers of m 2 /T 2 as we did in Section [| thereby obtaining additional contributions to D, 
B, and A. Doing so here, however, is not particularly enlightening so we will simply work 
with the expression given in Eq. ( |1.1| ). 

In analyzing the parameters in the tree-level potential of Eq. ( |2.1[) that lead to a strong 
first order EWPT, we perform numerical scans while implementing the various physical 
considerations discussed above. In doing so, we utilize the full finite temperature effective 
potential of Eq. (2.27) without resorting to an expansion in m 2 /T 2 . In order to interpret 



the results, however, it is useful to carry out an analytic study, relying on approximations 
justified by our full numerical studies. To that end we first note that our numerical study 
indicates that the effects of the BT 2 ip term in Eq. ( |4.1| ) are relatively unimportant compared 
to the tree-level cubic term proportional to e. Although both B and e depend linearly on 
the parameters a\ and 63, the factor of sin a/12 in B suppresses their effect, especially when 
the value of sin a is small near the critical temperature, as we find is typically the case. 



Thus, dropping the linear term in Eq. (44-) is a well-justified and helpful approximation. 

We now derive conditions on the parameters such that the condition for a strong first 
order EWPT (following from Eq. fll.2j ) and making use of the projection of v(T c ) over the 
SU(2) L direction g |l|, [D| ||) 



V ^-V2cosa c (^)>l (4.3) 



obtains. The values of ip c , a c , and T c are determined from the two minimization conditions 

dV cS ((p min ,a min ,T) 



dip 

dVcsifmrn, Q-min, T) 

da 



(4.4) 
(4.5) 



and the condition that defines T c , namely, that the potential at the minima of broken and 
unbroken electroweak symmetry as the same value: 

V c s{iPmin, Omin / 7r/2, T c ) = V e g((f min , 7r/2, T c ) = V , (4.6) 

where Vo is the high-T minimum defined in Section [||. 

This procedure is illustrated in Fig. [I] and |2[ The model we consider in Fig. [I] does 
not develop a singlet vev prior to the EWPT, and is characterized by the parameter set 
(ai, a-2,-1 63, &4, Xo, A) = (—933,0.69,356,0.53,55,0.4), with the parameters having mass di- 
mension 1 in units of GeV. The model we study in Fig. [2] features instead a non-zero singlet 
vev at T > T c , and is defined by (ax, a 2 , b 3 , 64, x , A) = (-88,0.765,-100,0.46,300,0.3). 
In the upper left panels we show the value of the potential as a function of the angle a 
at the location of the minimum in the radial ip direction, defined by Eq. (|4.4|). In the 
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Figure 1: Upper panels: The potential V e s (if m i n (a), a, T) (left) and the location of the minima 
in the (v/^/2,x) plane (right) for a model featuring the origin as its high-T vacuum, for various 
values of the temperature T. Lower panels: A contour map showing iso-level curves of constant 
values of the potential at T — 65 GeV (left) and at T = (right). The red arrows point towards 
smaller values of the potential. The black dots indicates the location of the absolute minimum of 
the potential. 



upper right panels, we indicate the location of the corresponding valleys of minima on the 
(v/y2,x) plane. The two lower panels show iso-level contour maps illustrating the behav- 
ior of the potential in the (v/y/2,x) plane at high temperature and at T = 0. The red 
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v/sqrt(2) [GeV] 



Figure 2: Upper panels: The potential V e s((p m i n (a),a,T) (left) and the location of the minima 
in the (v/y2,x) plane (right) for a model featuring a non-zero singlet vev at T > T c , for various 
values of the temperature T. Lower panels: A contour map showing iso-level curves of constant 
values of the potential at T = 340 GeV (left) and at T = (right). The red arrows point towards 
smaller values of the potential. The black dots indicates the location of the absolute minimum of 
the potential. 



arrows point towards smaller values of the potential. The black dots indicates the location 
of the absolute minima. 

Let us first consider Fig. [l|. From the upper left panel we see that when T = 65 GeV, 



for a/n < 0.062 the potential does not develop a minimum in the radial direction, hence the 
discontinuity in the resulting curve. The absolute minimum is at the origin, as indicated by 
the black dot in the right panel. Following again the curve at T = 65 GeV, we see that for 
0.062 < a/ir < 0.13 there exist minima in the radial direction, whose location is shown by 
the double-dotted-dashed line in the right panel. However, the values of the potential along 
that line are always larger than the value at the origin, as indicated by the same line in 
the left panel. Therefore, the absolute minimum is at the origin (the black dot in the right 
panel). Lowering the temperature, the potential on the line of radial minima decreases, 
until, at T = T c (dot-dashed line) its value at a = a c ~ 0.093-7T equals the value at the 
origin. In the (v/y/2,x) plane, upper right, the corresponding point at <p m i n (a c ) is shown 
by the black dot on the dot-dashed line. Finally, the black line, corresponding to T = 0, 
indicates that the minima in the radial direction exist now in the range < a/n < 0.13, 
and that the vacuum of the theory corresponds to ao — 0.95 and to the black dot on the 
solid black curve in the right panel on the (v/V2,x) plane. As explained above, the two 
lower panels help figuring out the shape of the potential at T = 65 GeV (left) and at T = 
(right); the black dots indicate the location of the absolute minimum of the potential. 

Similarly, in Fig. [2], the minimum of the potential at T = 137 GeV is at x ~ 385 GeV 
and v = 0, i.e. at a/ir = 0.5, as indicated in the left panel. At the critical temperature 
(dashed line in the left panel) the value of the potential at a. c jix ~ 0.34 equals that at 
a/ir ~ 0.5. The location of that minimum in the (v/y/2,x) plane is shown by the black 
dot on the solid black line in the right panel. For temperatures below T c the location of 
the minima in the radial direction doesn't change visibly on the (v/y2,x) plane, but the 
value of the potential at the true vacuum keeps decreasing, as indicated by the solid black 
line in the left panel. 

We now solve Eqs. (|4.4Q and (p~q) analytically and study the resulting potential 
V e g[(p m i n (a),a,T c (a)] as a function of a as above to find the minima of the potential. 
It is useful to consider first the case for Vq = 0. The solution to Eqs. ( |4.4j ) and ( [4.6D then 
yields 

where the "0" subscript denotes case for Vo = 0. As indicated by Eqs. (^4.2[), the values 
of e, e, and A are explicit functions of a while T c depends on a implicitly. The value of 
a = a c ^ 7r/2 that minimizes the potential is obtained from a numerical examination of 
Vctt [v?min(a), a, T c (a)] and determines the value of the ratio in Eq. (|4.7|), as described in the 
discussion of Fig. |]. In order to prevent washout of the baryon asymmetry by electroweak 
sphalerons, the projection of (f c /T c along the SU(2)l direction must be > 1 



V2cosa c (^)>l (4.8) 



leading to the result anticipated in Eq. ( |L4 

We analyze the case when Vo ^ by recalling from the analysis of Section || that 
one has [Vol/T^ << 1 for phenomenologically realistic parameters in the singlet potential. 
Thus, we solve Eqs. (|4.4j) and ( |4.6j ) by expanding in powers of | Vjj | / about the Vq = 
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solution and obtain, to leading non-trivial order in |Vo|/r^, 

+ ••• , (4.9) 



1+7^ 

c . 



where the + • • • indicate higher order terms in the | Vq \ /T 4 expansion and where 

= [3e-2~X( i p c /T c ) ][e-e/T c }-16~XD /^A ~ 4 
7 2e(e - e/T c ) - 16AL» V^Jo 



Here, the value T c appearing in Eqs. ( O| , 4.10| ) is the value obtained for Vq = 0. From 



Eqs. (p~8l) and (^9|) we obtain the result in Eq. (fT~4j) . 

It is relatively straightforward to demonstrate that the coefficient 7 is positive for 
realistic parameter choices 8 , so that the presence of a non-trivial high-T minimum always 
strengthens the EWPT. We have explicitly verified and confirmed this analytic result with 
our full numerical study. Thus, for practical purposes, it suffices to study Eq. (1.4) for 



Vq = in order to determine the potential parameters for which one obtains a sufficiently 
strong first order EWPT. Before doing so numerically, it is instructive to express this 
condition in terms of the parameters in the potential: 

4£ , 5 M cos 4 a c - v^cosacsinac (ai cos 2 a c /2 + 63 sin 2 a c /3) T c > ^ (4 11) 

2Xq(T c ) cos 4 a c + a2 cos 2 a c sin 2 a c + 64 sin a c /4 

where we have neglected the non-analytic dependence on the background field as before. 
We now make several observations about Eq. (f4.1l| ). 



(i) In the limit of pure SU(2)^ scalar sector, one has a c = and Eq. (|4.11 ) reproduces 
the SM result of Eq. ([D| with the substitution A (T C ) « A (0) = m 2 h /2v^. 



(ii) The introduction of a singlet extension of the scalar sector can lead to an increase 
in the strength of the EWPT in one of four ways 9 : (1) the presence of the terms 
in the numerator of Eq. fl4.11| ) associated with the breaking of Z2 symmetry, so 
long as cos a c sina c (ai cos 2 a c /2 + 63) < 0; (2) a reduction in the magnitude of the 
denominator due to a negative value for a>2 (but still satisfying the condition that 
the potential is bounded from below); (3) a reduction of the denominator due to 



singlet contributions to the one-loop effective potential that renormalize Aq [ 23 1 ; and 



(4) singlet contributions to the parameter e that we have not shown explicitly. 

(hi) Singlet contributions to e occur under specialized choices of the singlet parameters 
that lead to a cancellation of the mass functions II as in the MSSM case with the light 
scalar top, leading to an additional ip^T term in the potential beyond that of the SM. 



This special case has been studied in Ref. [21] for models having Z2 symmetry, and 
it was found that introduction of a large number of singlet scalars [~ 10] is needed to 
bring about a sufficiently strong first order EWPT for realistic masses for the SM-like 
Higgs boson. 



8 For instance, in the SM limit 16AD > e 2 . 



9 As noted in the introduction, choosing the parameters to give Vb < leads to a fifth avenue for 
strengthening the EWPT; the conditions for achieving this situation have been discussed in Section [| 
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High-T Trivial Singlet v.e.v. Non-zero High-T Singlet v.e.v. 




Figure 3: The correlation between the order parameter v c /T c and the critical temperature angle a c 
for models featuring a strongly first order EWPT for the parameter space scan outlined in Tab. |^. 
Black square points correspond to models also consistent with Higgs searches at LEP 2. The left 
panel shows models with a vanishing singlet vcv at T > T c , while that to the right models with a 
non-zero singlet vev at T > T c . 



(iv) The effect of singlet renormalization of A as studied in Ref. [23] is suppressed by a 
factor of 1/ 647T 2 , so a large value of the parameter 02 that enters this renormalization 
is needed to bring about an appreciable effect. Since the value of ai contributes to 
the degree of mixing between the singlet-like and the SM-like mass eigenstates and 
also governs the magnitude of T{hi — ► hjhj), a value of |<X2 1 that is large enough 
to produce substantial renormalization of Ao will have noticeable phenomenological 
consequences (see Sec. ^). The most dramatic effect, when nil > 2m,2, is a large 
suppression of the branching ratio for the decay of the lightest scalar to ordinary 
SM decay modes, as the largest decay width will correspond to a pair of singlet-like 
scalars. 

(v) In contrast, mechanisms (1) and (2) require no special tuning of parameters to yield 
a sufficiently strong first order EWPT when only one light scalar singlet is present. 
Not surprisingly, most of the models we find in our numerical scan correspond to 
these two cases. 

Before proceeding to low-energy probes of the parameters relevant to the EWPT, we 



consider Eq. (4.11) in the limit of small sina c . As indicated in Fig. ||, a strong first order 
EWPT occurs quite readily for small sina c , especially in the case where the singlet acquires 
a vev before the EWPT (right panel). Even in the left panel we find that a large fraction 
of points have cosa c ~ 1, although this is hardly visible and quantifiable from the figure, 
where the points are all aligned in the upper part of the plot. 
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Figure 4: As in Fig. g, for x and A (T = 0). 



Retaining only terms to linear order in sina c we obtain from Eq. ( [4.11 ) 

2Esm 1 a\ tana 



Ao 



2V2 T c A 



>1 



(4.12) 



Thus, one sees that in this regime, a relatively large negative value of a\/T c is needed for 



substantial enhancements of the EWPT compared to the SM situation. From Eqs. ( 2.1C - 
2.12j ), a\ also governs the degree of mixing between the singlet and SU(2)x neutral scalar 
in the case of small 02- Thus, for sufficiently light singlets, Higgsstrahhmg studies at e + e~ 
colliders can probe the parameter space needed for a strong first order EWPT. To see this 
situation explicitly, we consider the regime in which << ^ and express the LHS of 
Eq. (4.12) in terms of the mixing angle 9 and neutral SU(2)^ scalar mass: 



AEsm 



4 



2 tan 9 tan a c 
1 - tan 2 9 



> 1 



(4.13) 



For light /i 2 , large values of tan# are excluded by LEP searches for a light scalar with 
reduced couplings to the Z° boson. Consequently, we observe in Fig. |3| that imposing 
the LEP constraints significantly reduces the number of parameter space points that lead 
to a strong first order EWPT. In Section ||, we explore the implications of present and 
prospective Higgsstrahlung studies without resorting to the small sina c and n\ » [i 2 s 



limit. Nevertheless, the limiting case given in Eqs. ( 4.12 , 4.13 ) provides useful guidance for 
interpreting the more general numerical studies. 

It is also instructive to analyze the dependence of the LHS of the inequality ( [4.8| ) on the 
parameters in the tree-level potential. In Fig. ^, we show the values of the zero temperature 
singlet vev xq and SU(2)i quartic coupling Aq(0) that lead to a strong first order EWPT. 
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Figure 5: As in Fig. [3, for the mixing parameters a\ and 02. 



The points in Fig. [| were generated by scanning over the potential parameters as discussed 
above, and requiring a strongly first order EWPT; points in black are consistent with the 
LEP Higgs search constraints on the (?7ifc,£?) plane (A; = 1,2), while the black and light 
colored points together indicate the results of our scan over all models. The results are 
shown for the case when both the high-T singlet vev is zero (left panel) or non-vanishing 
with Vq < (right panel). We observe that models with small values of Ao(T = 0) are 
excluded by LEP since they lead to exceedingly small values for the mass of the SM-like 



neutral scalar [see Eq. ( 2.10 )1. Moreover, models having a substantial zero-temperature 
singlet vev and a strong first order EWPT occur quite abundantly, even after the LEP 
constraints are implemented. We also notice that, unlike the case where the singlet acquires 
a vev before the EWPT, when the EWPT proceeds from zero vev for the singlet (left panel) 
various models feature very small values for xq and are LEP-compatible (the vertical string 
of points in the left part of the plot). For those models, fj^ is typically very large due to 
the ciiVq/(4xo) term in Eq. 2.10| , and the mixing between the two scalars is suppressed 

In Figure || we show the distribution of parameters a\ and 02 that couple the SU(2) l 
and singlet Higgs and that lead to a strong first order EWPT. We note that in nearly 
all cases, the Z2 symmetry-breaking parameter a\ is negative, as one expects from the 
discussion of Eq. (4.11). Models with 02 of either sign can lead to a strong first order 
EWPT whether or not the LEP constraints are implemented. Note that when \a\\ is large, 
a positive value for a 2 is more likely to be consistent with the LEP constraints. This trend 
can be understood qualitatively from the dependence of on these two parameters; to 
avoid large mixing between the singlet and SU(2)^ neutral scalars, the value of 02 must 
be relatively large and positive to cancel a large negative contribution to /x? from a large 
negative a\ when xq differs significantly from zero. When both a\ and 02 are negative, 
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Figure 6: The (A, 02) and (x, 02) planes for models with a\ = 63 = 0. 



however, the constraints on large mixing can be relaxed by making fi 2 s large by suitable 
values of 63,4, a situation corresponding to the points in the lower left quadrants of both 
panels in Fig. B We also find that in models with Z2 symmetry, a strong first order EWPT 
is generally consistent with the LEP constraints only when 02 < 0, as indicated in Fig. |6|. 
Again, this behavior is consistent with Eq. (4.11) since a negative value of 02 can reduce 
the denominator of LHS without engendering large singlet-doublet mixing. 

Finally, we show in Fig. the dependence of \[2 cos a c <p c /T c = v c /T c versus T c . Here, 
we see that imposition of the LEP constraints has a substantial impact. In the absence 
of such constraints, the value of T c leading to a strong first order EWPT can vary from 
as much as ~ 200 GeV to 10 GeV, whereas imposing the LEP bounds favors critical 
temperatures of order 100 GeV. Hence, for the Vb / case, our expansion in | Vb|/T^ is 
good for the LEP-allowed models. 



5. The zero-temperature phenomenology 

As the discussion of Section |^ indicates, searches for the Higgs boson and studies of its 
properties can provide important information about the viability of a strong first order 
EWPT in the presence of an extended scalar sector. Here, we discuss the present implica- 
tions of EWPO as well as the possible tests of this scenario that might be carried out at 
the LHC and at a Linear Collider. 

Electroweak Precision Observables & Singlet Scalars 

We first consider the impact of EWPO on singlet models that also lead to a strong first 
order EWPT. We address this issue by computing the scalar contributions to the gauge 
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Figure 7: As in Fig. pi for the order parameter versus the critical temperature. 
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Figure 8: Feynman diagrams for scalar one-loop corrections to the gauge boson propagators. 

boson propagators (Fig. |8|) in the Feynman-'t Hooft gauge and compare the corresponding 
effects with those from the SM Higgs on EWPO. Although the masses are generally 
of order the electroweak scale, it is possible to characterize their effects on the Z°-pole 
observables and VF-boson mass with the oblique parameters S, T, U, and V. The effects on 
lower-energy precision observables, such as the weak charge of the cesium nucleus obtained 
from atomic parity violation [^] or the parity-violating asymmetry in M0ller scattering 
1 33] require that one take into account the g 2 -dependence of the gauge boson self energies 
that is often approximated using the additional oblique parameters W, X, etc.. In the 
present analysis, we will neglect these higher derivative terms, saving a more complete 



analysis to a future publication p4|| . Note that the introduction of a real neutral scalar 
does not affect the propagator functions 1T 77 and IIz 7 , so it suffices to consider only the 
effects on IIww and ILzz in the present situation. 

A simple indication of the difference between the effects of the SM Higgs and those 
of the mixed states h%2 on the precision observables can be seen from the expressions for 
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scalar contributions to the T parameter in the doublet-singlet scenario as compared to the 
SM case: 



T — Tsm 




(5.1) 



where s 2 = sin 2 6 'w(Mz) gives the weak mixing angle in the MS scheme at the t'Hooft 
scale \i = Mz, c 2 = M^/M^, and rrih is the reference value of the Higgs boson mass in 
the SM. Note that the effect of the singlet enters both through the presence of the mixing 
angles in the prefactors as well as the values of the masses m the expressions that they 
multiply. The expressions for S — Ssm and U — Usm are analogous, but too cumbersome 
to reproduce here since they involve evaluation of the loop integral functions 



F n (ml,ml, 



q 2 ) = I dx x n In [(1 — x)m 2 + mij — x(l 
Jo 



(5.2) 



In order to identify the EWPT-viable models that are favored by EWPO, we have 
performed a fit to the precision observables identified in the Review of Particle Physics 
(RPP) [|]] for a fixed value of = 114.4 GeV and obtained the values of AO = O — Osm, 
O = S, T, U. To carry out this fit, we have employed the GAPP program |}5| that has 
been used in arriving at the best-fit values in the RPP, but using the value of the top quark 
mass rat = 170.9 ± 1.8 given in Ref. [pq]. The results give the best fit values and standard 



errors 



10 



T -T S m = -0.111 ±0.109 
S-Ssm = -0.126 ±0.096 
U -U SM = 0.164 ±0.115 



(5.3) 



We define a Ax 2 via 



A X 2 = £(AO ? , - AO?) (AOj - AO®) , 



(5.4) 



where AO? denotes the best fit value for the differences in a given oblique parameter from 



pij 



its SM reference value; afj = aiPijdj, with dj denoting the errors in Eq. (|5.3|) ; and the 
correlation matrix is 

/ 1 0.866 -0.588 \ 

0.866 1 -0.392 . (5.5) 
y -0.588 -0.392 1 J 



10 The GAPP routine includes SM Higgs radiative corrections explicitly. Consequently, the SM Higgs 
contributions to the oblique parameters for a fixed value of nih must be subtracted from the singlet- doublet 
scalar contributions when interpreting the results of the fit. 
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Figure 9: AO, for O = T,S, U, and A% 2 as a function of the SU(2)-like Higgs mass m\ for 
three values of the singlet- like Higgs mass 1112 and the mixing angle sin#, respectively (771 2 = 
50 GeV,sin# = 0.5), black solid line, (777,2 = 50 GeV,sin6> = 0.1), red dashed line, and (7772 = 
500 GeV, sin# = 0.5), green dot-dashed line. 



The 95% C.L. ellipsoid in the space of AOj corresponds to taking A% 2 < 7.815. The 
Ax 2 for the SM case of a single Higgs scalar with mass given by the direct search lower 
bound (mh = 114.4 GeV) is obtained by setting the AO{ = 0, yielding A% 2 = 3.546. We 
consider models with the singlet scalar extension to be consistent with the EWPO if they 
generate values for the AOi lying within the 95% C.L. ellipsoid. The contributions from 
the extended scalar sector to S, T, and U - and, thereby, to electroweak observables - is 
governed by three parameters: the two masses mi 2 and mixing angle 9. As indicated by 
Eq. (5.1), the dependence of the the oblique parameters on the masses an mixing angle is 
non-linear, so we would not expect a three-parameter fit to EWPO using S, T, and U to 
be identical to a three-parameter fit using mi^ (or In 777,^2) and 9. Nevertheless, we expect 
the 95% C.L. ellipsoid for the oblique parameters to provide a reasonable indication of the 
consistency between EWPO and the extended scalar sector models. 

We show in Fig. |9| our results for the differences AO between the oblique parameters 
T, S and U in the present model and their reference SM values, as well as the quantity 
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Ax 2 - We consider the variations of all of these quantities as a function of the SU(2)-like 
scalar mass, mi, for three different combinations of the singlet-like mass rri2 and of the 
mixing angle sva.9. In particular, we consider two cases with large sin# = 0.5, the first one 
with small ni2 = 50 GeV 11 and the second one with large 7772 = 500 GeV, and a third 
model with 7772 again set to 50 GeV, but smaller mixing angle, sinf? = 0.1. 

It is interesting to compare the effects of AS and AT - and their dependence on mi,2 
and 9 - on the overall A% 2 (the effect of AU is small 12 ). In the case of AT, each of the last 
two terms in Eq. (5.1) gives a negative contribution. Moreover, the magnitude of each term 



grows as the scalar mass is increased. Thus, AT tends to favor larger scalar masses and 
smaller mixing when 777-2 < mi- In contrast, AS favors lighter scalar masses and, therefore, 
larger mixing in the 777-2 < m\ regime. In the competition between the two effects, AS 
tends to dominate, driving the Ax 2 toward lighter scalars with larger mixing. 

It is already known that EWPO favor a light SM Higgs, with the minimum x 2 occurring 
for m/j below the direct search bound of 114.4 GeV. In the case of the doublet-singlet 
scalar sector considered here, one would expect the precision data to favor the presence of 
a second mass eigenstate that is light, and that has relatively large couplings to the SM 
gauge bosons. This trend, indeed, emerges from our numerical analysis, as indicated in 



Fig. 10, where we show the contours of constant Ax 2 in the (mi, 777-2) plane for two different 
values of sin#. For small sin# (left panel), the effect of the singlet-like scalar is weak, with 
the Ax 2 being governed almost entirely by the mass of the SU(2)-like scalar. Including 
a very light singlet-like scalar improves the A% 2 , but not substantially. In contrast, for 
large sin#, corresponding to mass eigenstates that are nearly equal mixtures of singlet and 
SU(2)/, components, one can improve the Ax 2 from the SM reference value if at least 
one of the two scalars is light. The LEP Higgs searches still exclude much of this region, 
since the mixing is nearly maximal. Nevertheless, the EWPO allow for both scalars to be 
considerably heavier than 114.4 GeV (up to about 220 GeV) in the presence of maximal 
mixing. 

From these considerations, we expect that the impact of EWPO on models featuring 
a strongly first order phase transition will be to limit mi to values smaller than around 
200 GeV; large values of the mass of the singlet-like scalar are EWPO-viable as long as 
the mixing with the SU(2)-like scalar is suppressed. We confirm this trend comparing 



the upper and lower panels of Fig. 11, where we show the same models we showed in the 
previous plots in the (7771,7772) plane, without (upper panels) and with (lower panels) the 
EWPO 95% C.L. constraint. Imposing the latter eliminates nearly all models with a very 
light singlet-like scalar and heavy (mi > 200 GeV) SU(2)-like scalar. 

Collider Implications 

As a starting point of our analysis of collider implications, we refer to Eqs. ( 2.10| - 2~T5| ) 



11 We employ here this value of m-2, that would be excluded by the LEP 2 Higgs search results, for 
illustrative purposes only, to show the maximal effects on EWPO of having a light singlet Higgs. 

12 The U parameter is defined as the difference of finite differences involving ffn and II33 [^. In the 
limit of vanishing weak mixing angle, the neutral scalar contributions to U vanish, so it is not surprising 
that the scalar contributions to this parameter are quite small. 
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that indicate how the masses and mixing of the fields characterizing fluctuations about the 
zero temperature vevs depend on the parameters in the full potential. We first consider 
the possibility that Higgs production and decays at the LHC could be modified by the 



singlet extension as recently discussed in Ref. [20] and elsewhere j38|]. If mi > 2m2, then 
the decay hi — > ^2^2 becomes kinematically allowed, and the branching ratio (BR) of hi 
to the conventional SM final states will be reduced. The presence of this additional Higgs 
decay channel can provide interesting signals as four fr-jets or 2r + 26; this possibility has 
been explored in various supersymmetric extensions [39|. Conversely, if m<i > 2mi, then 
double SM-like Higgs production becomes possible if hi is first produced, leading to similar 
exotic final states. However, the production rate of /12 is suppressed relative to hi by the 
mixing factor tan 2 6 < 1. 

Fig. 11 illustrates the distribution of points in the (7711,7712) plane that give rise to a 
strong first order EWPT, with (lower panels) and without (upper panels) imposing the 
EWPO constraint. As before, the distribution of all models is given by the light and dark 
points, while those consistent with LEP searches are only those in black. Notice that we 
find several models featuring a large value of the singlet-like Higgs mass eigenstate m-z' for 
those models, since 7712 3> T c , finite temperature corrections play a negligible role in making 
the EWPT strongly first order. Rather, the driving effect in the Higgs effective potential is 
played by tree level terms involving the field S, as discussed in Sec. |||. Points lying below 
the 7772 = mi/2 blue line correspond to kinematically allowed hi — > /i2^2 decays, while 
those above the blue 7772 = 2mi line indicate models in which the h% — ► hihi channel is 
open. 

If one does not impose the EWPO constraints, a substantial fraction of the EWPT- 
viable models also give rise to either of these two open channels, indicating the possibility 
that such models may be probed with Higgs studies at the LHC. A detailed analysis of such 
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Figure 11: As in Fig. ||, for the physical Higgs masses (mi corresponds to the mostly SU(2)-like 
mass eigenstate, and 7712 to the mostly singlet-like mass eigenstate). The two upper panels show 
all models, while the lower panels show the models that are not ruled out at 95% C.L. by EWPO. 
For models in the lower- right part of the plot (below the lower blue line), the decay hi — > /12 h,2 is 
open, while for models above the upper blue line the decay hi — > hi hi is open. 



searches will appear in a companion paper [p0|| . It is striking, however, that imposing the 
EWPO constraints closes - for a large fraction of the models in our scan - the window of 
hi — ► /12/12 decays for models also compatible with LEP 13 . In contrast, the decay h^ — > h\h\ 
is generically open for a large set of models, provided that m\ < 160 GeV. 

As a special case, we focus on models having a Z2 symmetry before SSB, given by the 

13 Models in this region still occur, but they do not appear in the plot due to the resolution of our scan. 
Their existence is more apparent in the special case of Z2 symmetry discussed below. 
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Figure 12: Right: The (mi, 7712) plane for models featuring a\ — 63 = 0. Left: for the same 
models, the ratio of the decay width of the SU(2)-likc Higgs into two singlet, over the total decay 
width of the SU(2)-like Higgs into Standard Model particles. The two upper panels show all models, 
while the lower panels show the models that are not ruled out at 95% C.L. by EWPO. 



condition a\ = 63 = 0. Such models are particularly interesting since they also admit the 
possibility that h,2 is the dark matter particle, since it cannot decay to any SM particles 
when xo = and there is no mixing with the SM-like Higgs [25]. The corresponding 
distribution of EWPT-viable models in the (mi, 1712) plane is given in the left panels of 
Fig. [l^. Notice that here we consider a special slice of the parameter space, as we scan 
only over four of the six parameters in Tab [2|, setting a\ = 63 = 0. Thus, the ensemble of 
~ 10 4 models we find differs from those we considered above. 
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An important feature emerging from this scan is that the LEP-consistent models always 
lie in the h\ — * h%h% allowed region. Even after EWPO constraints are imposed, a large set 
of models exist for which this decay channel is open. For these models, one would expect a 
reduction in the BR of hi to SM modes. The degree of this reduction depends on the rate 

#122 



where 



T(h 1 ^h 2 h 2 ) = f^\ll-Aml/m\ , (5.6) 



5122 = — ^— {a 2 {c e - 2s e ) + Q\s e ) + (s% - c g s e ) + 

+(6 3 + 36 4 a;o) c 2 e s e , [sg = sin#, c e = cos 8}. (5.7) 

We note in passing that if T(hi — > h 2 h 2 ) > y sm , models with a light SU(2)-\ike Higgs 
featuring mi < 114.4 GeV can be compatible with the LEP exclusion limits [jTJ] [see 
Eq. ( |2.17| )1. The rate for decay to SM particles T(hi — > SM modes) depends on the final 
states that dominate for a given value of mi. For light hi, its SM decays are dominated 
by the bb final state. In Ref. |20f| , the BR reduction was studied for this regime under the 
implicit assumption that xq = (sin# = 0). In this case, fj, 2 = b 2 /2 + a 2 v 2 /2 with b 2 an 
independent parameter, so that increasing the coupling a 2 for fixed b 2 affects the rate in 
Eq. ( |5.6| ) through both the overall coupling strength as well as the mass of the h 2 . 

Here we consider the more general situation (compared to the foregoing) in which 
xo 7^ is an independent parameter. In the right panels of Fig. [l^ we show the ratio 
T(hi — > h 2 h 2 )/T(hi — * SM modes) for the EWPT-viable models, before (upper panel) and 
after (lower panel) imposing EWPO constraints. For light hi (114 < mi < 180 GeV), this 
ratio can be quite large, leading to a suppression of the SM BR by an order of magnitude 
or more. For heavier hi, the decay competes with decays to massive gauge bosons, leading 
to a reduction in the SM BR from one to more than ten percent, while for mi > 300 GeV, 
the reduction is generally less than ten percent and for many models can be well below one 
percent. However, as can be seen in the lower right panel, these models are incompatible 
with EWPO constraints. For models consistent with EWPO at the 95% C.L. we always 
expect a BR suppression at at least the percent level. In short, if a SM-like Higgs with 
mi < 200 GeV is discovered at the LHC, one would expect an observable reduction in 
its SM BR in the presence of a light singlet-like scalar that gives rise to a strong first 
order EWPT. For those models giving a BR reduction of < 10%, Higgs studies at a Linear 
Collider would likely be needed to observe this suppression. Conversely, the absence of a 
BR reduction for heavier SM-like Higgses would point to models that do not have a Z 2 
symmetry before SSB or have a heavier singlet-like scalar. In this case, one would look for 
signatures either through exotic final states at the LHC, as discussed above, or through 
future Higgsstrahlung studies at a linear collider. 

In scenarios having no Z 2 symmetry, there exists a substantial fraction of EWPT- 



viable models that lie between the two blue lines in Fig. 11 and that imply no change in 
the Higgs decay modes at the LHC. For these models, it is interesting to ask whether future 
Higgsstrahlung studies at a Linear Collider could provide a useful probe. To that end, we 
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Figure 13: The (mi, ma) plane for all models considered in the various scans, in the region where 
the decays — > ^2.1^2,1 are kinematically forbidden. In the right panel we removed models that 
did not satisfy the EWPO constraints at 95% C.L. 



show in Fig. 13 the distribution of points where h\2 — * /12, 1^2,1 decays are kinematically 
forbidden, focusing on those having a small value for £ 2 that characterizes the reduction in 
the Higgsstrahlung rate. Points with m<i < 114.4 GeV that are not ruled out by the LEP 
Higgs search generally have £ 2 << 1. For mi 2 > 114.4 GeV, the overwhelming majority 
of models have £ 2 > 0.01 and generally correspond to large values of a\ as implied by 
Eqs. (|2.1C|J2.12D . Those points are shown with orange dots in the figure. Points with 
0.001 < i 2 < 0.01 (dark green diamonds) and £ 2 < 0.001 (black squares) indicate the 
Higgsstrahlung sensitivity needed to probe these models. To the extent that Higgsstrahlung 
studies at a future e + e~ collider can search for neutral scalars in this mass range with 
Higgsstrahlung rates up to 100 times smaller than for a SM-like Higgs, one could probe 
a large class of EWPT-viable singlet models that would not be readily accessible at the 
LHC. Finally, the panel to the right illustrates the impact of imposing EWPO constraints 
on the models at hand. As elsewhere, including these constraints considerably reduces the 
space of viable models. 



6. Conclusions 

Searches for the Higgs boson at the LHC and studies of its properties both there and at a 
future Linear Collider will be an important thrust at these colliders. In this study, we have 
analyzed the implications for Higgs boson phenomenology and the EWPT of the simplest 
extension of the Standard Model scalar sector containing one additional gauge singlet scalar 
field. Our goal has been to identify general features of more complicated models with 
gauge singlet scalars, such as extensions of the MSSM in which the ^-term is generated 
dynamically as the vev of such a scalar field. Although some details of Higgs phenomenology 
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as it relates to the EWPT will depend on specific model realizations, we believe we have 
outlined general trends that one would expect in most scenarios of this type. In doing so, 
we derived a simple criterion that allows one to predict, starting from the parameters in 
the scalar potential, whether or not the singlet field acquires a vev prior to the EWPT. We 
also obtained a compact expression [Eq. ( |1.4| )] for the condition of obtaining a strong first 
order EWPT as needed to prevent washout of the BAU by electroweak sphalerons, under 
approximations that were systematically cross-checked through our numerical analysis. 
The features that emerge from our study include: 

- We find that a strongly first order EWPT can occur quite readily in models where 
the singlet scalar obtains a non-zero vev (xq ^ 0). Choosing xq = (see e.g. [plf ) is 
a necessary condition for models that have a scalar dark matter candidate and forces 
the potential to have a Z2 symmetry. Conversely, models with a Z2 symmetry before 
SSB do not necessarily have xq = 0. 

- Models with singlet scalars can readily enhance the strength of the EWPT to allow 
viable electroweak baryogenesis through one of three tree-level effects: the presence 
of a sufficiently large negative coupling between (H'H) and S, thereby enhancing the 
numerator of Eq. ( |1.4j ); the presence of a sufficiently large negative coupling between 
(H^H) and S 2 , thereby reducing the denominator of Eq. ( |1.4[ ); or the choice of pa- 
rameters in the singlet scalar potential that generates a non-trivial high temperature 
minimum (before electroweak symmetry is broken). 

- The presence of a Z2-violating (if^i7)S* coupling leads to mixing between the singlet 
and neutral SU(2)i scalars, as does the ^-invariant (H^H)S 2 interaction when xq 7^ 
0. Higgs boson searches at LEP exclude models in which one or the other of these 
couplings is large enough to induce a strong first order EWPT when one of the scalar 
masses is well below the direct search lower bound. Future Higgsstrahlung studies 
at a Linear Collider could provide an additional probe of these effects for heavier 
scalars 14 . 

- Considerations of EWPO imply that EWPT-viable models in which the singlet-like 
scalar can decay to two SU(2)-like scalars occur more copiously than those in which 
the SU(2)-like scalar decays to a pair of singlet-like scalars. In principle, the models in 
which the singlet-like scalar decays to a pair of SU(2)-like scalars could be identified 
through the presence of exotic final states at the LHC, such as four 6-jets, 66r + r~, or 
6677. EWPT-viable models in which these decays are kinematically forbidden could 
still be tested with Higgsstrahlung studies at a Linear Collider. 

- In the special case of models having Z2 symmetry before SSB, all EWPT-viable 
models that are consistent with EWPO lie in the region where the SU(2)-like scalar 
can decay to pairs of the singlet-like scalars. For these models, one would expect 
to see large suppressions of the branching ratio to SM Higgs decay modes, as the 
branching ratio to the pair of singlet-like scalars would be considerable. 

14 In principle, studies of Higgs production via WW and ZZ fusion could also be used for this purpose. 
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In short, the possibility of testing the viability of producing the strong first order 
EWPT needed for successful electroweak baryogenesis in models with a small number of 
singlet scalars appears to be feasible using a combination of Higgs studies at the LHC and 
a future Linear Collider. Together with the next generation of "table top" searches for 
the electric dipole moments of the electron, neutron, neutral atoms, and nuclei that will 
provide powerful probes of new electroweak CP-violation, up-coming collider Higgs studies 
could help unravel the origin of baryonic matter. 
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